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Abstract - This paper presents a new type of an adaptive fuzzy operator for detection of 
isolated abnormalities, and enhancement of raw ultrasonic images. Fuzzy sets used in 
decision rules are defined for each image based on empirical statistics of the color 
intensities. Examples of the method are also presented in the paper. 


I. Introduction 

Ultrasonic testing is one of the techniques used for nondestructive characterization of 
materials. In this work our interest is in one of its variations called contact scanning, where one 
front surface and two back surface ultrasonic pulses obtained using the pulse-echo configuration 
are digitized and stored at every scan location. These digitized pulses are used to compute 
ultrasonic properties versus frequency like phase velocity, cross-correlation velocity (not a 
function of frequency), reflection coefficient, attenuation coefficient and others [2-4]. Ultrasonic 
properties can be used to form ultrasonic images where values of a given property (at any 
frequency within transducer bandwidth) are represented by corresponding gray or color scale 
values. On a raw ultrasonic image each point (pixel) represents a single location for which 
ultrasonic data were collected. An example of a raw ultrasonic image representing a phase 
velocity 1 at 60 MHz is shown in Figure 1. Ultrasonic images from contact scans are typically 
enhanced using linear or cubic interpolation. Enhanced images are used to determine areas which 
correspond to undesired properties of a tested material. Figure 2 represents an enhanced ultrasonic 


1 Explanation of phase velocity is given in the Appendix. 
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image obtained directly from the raw image shown in Figure 1. Note three red and two blue 
"circles" in the image (dark circles on black and white illustrations). Each of these circles 
corresponds to a single point on the raw ultrasonic image; they are examples of "undesired 
phenomena" for further image analysis. In general, while constructing an ultrasonic image we are 
interested in global trends, thus properties of separated points are irrelevant. These "undesired" 
isolated points on the raw ultrasonic image should be identified and filtered before the image is 
enhanced. Figure 3 shows an enhanced image created from a filtered ultrasonic image of Figure 
1 . 

Existence of the undesired points on the ultrasonic images can be attributed to: 

1. Corruption of the ultrasonic pulse data. It usually occurs during digitization of 
analog data and has random nature. There are several types of these pulse data 
abnormalities, and only few of them can be detected by a direct analysis of digital 
pulses. 

2. Ultrasonic pulse data is collected and digitized correctly, but properties of the 
tested material at these isolated points are "significantly” different from its 
neighbors. 

3. Intrinsic uncertainties in data used for calculation of ultrasonic properties, like 
total difference in phase between first and second back reflection, used to calculate 
phase velocity. 

The aim of this paper is to construct adaptive fuzzy operators for automatic detection of 
isolated suspicious points on an ultrasonic image. These points once detected can be either 
removed to facilitate creation of smoother enhanced image, or further investigated to determine 
the nature of their origin. 

A method of fuzzy operator synthesis, presented in this paper, employs usage of statistical 
data obtained from a raw ultrasonic image, like histograms of pixel intensity or histograms of 
intensity differences between neighboring pixels. The data are processed and used to generate 
fuzzy sets representing desired versus undesired image contents. 


II. Structure of the Fuzzy Operator 

Let us start with a description of a fuzzy operator for image processing similar to the one 
presented in [1]. 


A. Generalized description of fuzzy variables 

Each pixel of an image is processed by fuzzy rules applied to a set of pixels belonging 
to a rectangular window (neighborhood) centered on that pixel. As the window scans each pixel 
of the image, one at a time, a new pixel of the resulting image is generated by means of fuzzy 
reasoning. Fuzzy variables can be defined by using any function (or relationship) of interest for 
the processing and can be expressed as followings: 


2 



(1) 


V 2 = 4^1 

rl n 1 X m 

o =fiP c J>') 


where v„ v 2 , ..., v n are the input variables, o is the output variable, /?,, p 2 , p m are luminance 
(intensity) values of the pixels in the window, p c and p\ indicate luminance of the pixel in the 
center of the window before and after processing, respectively. 


B. Fuzzy Rulebase 


The reasoning structure uses two types of rules: THEN -rules and ELSE-rules. Each 
THEN-rule includes n antecedents linked by fuzzy "AND" logical operators and only one 
consequent. The antecedents are related to input variables and the consequent to the output 
variable. The overall rulebase is composed of many THEN-rules and only one (global) ELSE- 
rule. The following rulebase contains r THEN-rules: 


IF ( v, is A„) AND ... AND ( v n is A al ) THEN (o is £,) 
IF ( V! is A n ) AND ... AND ( v n is A^) THEN io is B 2 ) 


IF ( Vj is A u ) AND ... AND ( v n is AJ 


THEN (o is B,) 
ELSE (o is C„) 


( 2 ) 


where A j} is a fuzzy set associated with variable v, in the y'-th rule, B } is the consequent set of the 
same rule and C 0 is the consequent of the ELSE-rule. The ELSE-rule assures that the inference 
is successfully executed even if no THEN-rules are fired. 


HI. Synthesis of the Fuzzy Operator 


A. Approach 

In the approach presented here we attempt to follow a human way of deciding which 
points on an a raw ultrasonic image might be considered "suspicious”. A typical decision is 
simple: 

A point is suspicious if its intensity is significantly different from the intensity of 
the surrounding points. 

We shall try to implement this "simple" decision in a fuzzy rulebase. The main difficulty is 
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establishing what "significantly different" constitutes. The difficulty stems from the fact that 
human decisions are more qualitative than quantitative, and are not easily translated into 
numerical values. 

We shall approach this problem by investigating statistical properties of image intensities. 
Figure 4 shows a histogram of pixel intensities of the raw ultrasonic image of Figure 1. Figure 
5 presents histogram of an absolute value of intensity difference between neighboring points. 
Underlying distributions, represented by these histograms, can be used to represent "suspicious" 
points. We can mark a point as suspicious if its intensity value is in a tail of the image intensity 
value histogram, and its minimum difference of intensity to neighboring points is in a tail of the 
image intensity difference histogram. Since distribution types (image intensity, intensity 
differences) are unknown and are dependent on an image content, standard methods of statistical 
analysis are not well applicable here. Thus, we shall investigate possibility of using these 
histograms to generate fuzzy sets representing points with "atypical" intensity values and fuzzy 
sets representing points in which intensities are "quite" different than their neighborhood. 


B. Creation of fuzzy sets 


Let us introduce a normalization operator, N: 


mx)] = 


Ax) -min. fa) 
max/(;c) -min/(jc) 


(3) 


After normalization, function f(x) takes on values from the interval [0,1]. Let s denote pixel 
intensity value, and d denote a minimum difference between the given point intensity and 
intensities of its neighboring points. 

The fuzzy set, A d , representing high-intensity-difference is created by taking the negation 
of the normalized inverted sum of an intensity value histogram, hjd)\ 

Ajd) = 1 -N[fh d (x)dx] (4) 

d 

Plot of the membership values of the fuzzy set A d versus intensity difference d is shown in Figure 

6 . 

The fuzzy set A s representing atypical-intensity is obtained by taking a minimum of the 
negated cumulative sum, and the negated inverted cumulative sum, of the image intensity value 
histogram, hjs). 

S 

H *{s) = f h(x) dx (5) 

J J S 
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( 6 ) 


Hjs) = J h(x) dx 

S 


A (s) = minil -N[H '( 5 )], 1 -N[H\s )] } (7) 

Plot of membership values of the fuzzy set A s versus intensity value 5 is shown in Figure 7. 

Using definitions of the above defined fuzzy sets we construct the following decision rule: 
The output of evaluating a rule is a certainty value 2 , in the 0 to 1 range, representing degree to 

IF (point’s intensity is a very atypical-intensity) 

AND ... 

(minimum difference of point's intensity is a very high-intensity-difference) (») 


THEN this point is suspicious 


which a point is "suspicious". For the purpose of filtering we can set some certainty level, so 
when it is met we can mark the point as the one to be filtered out. 


IV. Examples 

The algorithm was tested on phase velocity images for four circular ceramic samples, 
named SN2, SN3, SN4 and SN6 [5]. The first sample, SN2, was used for calibration, the other 
samples were used for testing. Figures 8 to 11 shows examples of modification of fuzzy sets by 
hedge operators. Two standard hedge operators, very and extremely , were tested. They are 
implemented by rising membership values of a fuzzy set to the power of 3 and 5 respectively. 
First two columns of Table 1 compare number of points selected as "undesired" (suspicious) 
using different certainty levels shown in first column, and high and extremely hedges shown in 
second and third column. In general, extremely is less permissive than very in selecting points, 
qualitatively however the points selected are usually the same. Samples SN3, SN4 and SN4 were 
evaluated using the hedge very. 


2 Certainty value should not be mistaken with probability. 
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Points selected 

Certainty value 

SN2 very 

SN2 

01 1 ^extremely 

SN3 

SN4 

SN6 

0.99 

1 

1 

1 

1 

1 

0.95 

5 

4 

5 

8 

7 

0.9 

8 

6 

7 

10 

9 

0.85 

11 

8 

16 

19 

12 

0.8 

14 

9 

19 

30 

17 

0.75 

17 

11 

24 

33 

22 


Table 1: Number of points selected for different samples. 


Figures 12 through 17 represent graphically how different confidence levels, using sample 
SN2, affect the quality of the filtered and then enhanced image. Based on these results level of 
confidence 0.9 was chosen. Figures 18 through 23 show comparison between raw enhanced 
images and images which were first filtered using our algorithm and then enhanced. The 
confidence level used for all of them was 0.9. Note changes in color scale: for filtered images 
the ranges are smaller. It is a highly desirable effect since it helps to better visualize 
microstructural gradients. 


V. Conclusions 

In this paper an adaptive fuzzy image processing operator was presented. The operator 
was primary developed for processing of raw ultrasonic images, however it can be easily applied 
in other image processing areas. The adaptiveness of the operator steams from the way the fuzzy 
sets used in the operator's fuzzy rulebase are generated. The sets are generated for each processed 
image individually using statistical data about image pixel intensities and intensities differences. 

The introduced, relatively simple adaptive fuzzy operator is quite robust in finding 
suspicious points in raw ultrasonic images. In the presented fuzzy operator detects only isolated 
points, it would be desirable, however, to make i f sensitive not only to isolated points but also 
to small areas. Other fuzzy rules can be added to the operator rulebase, for instance to improve 
treatment of border points. 

In presented examples the adaptive fuzzy operator was used for filtering of ultrasonic 
images. Other interesting are of its application is location of suspicious points which are than 
further processed by other means to investigate the origin of points abnormality. 
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Appendix - Ultrasonic Properties 


Notation 

t - time. 

/ - frequency. 

b,(t), b 2 (t ) - First and second back surface reflection. 

B,(f), B 2 (f) - Fourier transforms of first and second back surface reflections. 
X - sample thickness. 


Phase Velocity 

_ (2X)2 nf 

~ arg{B<f)) -arg(B 2 (f)) 


(9) 


where arg(B(f)) is an argument of Fourier transform of pulse B(t). 
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Figure 6: Fuzzy set: high-intensity- difference. 
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Figure 8: Fuzzy set: very high- intensity - 

difference. 
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Figure 7: Fuzzy set: unaverage-intensity-value. 
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Figure 9: Fuzzy set: very unaverage-intensity- 
value. 



Figure 10: Fuzzy set: extremely high-intensity- 
difference . 
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Figure 11: Fuzzy set: extremely unaverage- 

intensity- value. 






















